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ABSTRACT 
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Type IIP (plateau) supernovae are thought to come from stars with initial 
^ \ mass ~ 8—25 M that end their lives as red supergiants. The expected stellar end 

points can be found from evolutionary calculations and the corresponding mass 
^ | loss properties at this point can be estimated from typical values for Galactic 

| stars. The mass loss densities of observed supernovae can be estimated from ob- 

servations of the thermal X-ray and radio synchrotron emission that result from 
the interaction of the supernova with the surrounding wind. Type IIP supernovae 
are expected to have energy-conserving interaction during typical times of obser- 
\q \ vation. Because Type IIP supernovae have an extended period of high optical 

luminosity, Compton cooling can affect the radio emitting electrons, giving rise 
to a relatively flat radio light curve in the optically thin regime. Alternatively, a 
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high efficiency of magnetic field production results in synchrotron cooling of the 
radio emitting electrons. Both the X-ray and radio luminosities are sensitive to 
the mass loss and initial masses of the progenitor stars, although the turn-on of 
radio emission is probably the best estimator of circumstellar density. Both the 
mass loss density and the variation of density with stellar mass are consistent 
with expectations for the progenitor stars deduced from direct observations of 
recent supernovae. Current observations are consistent with mass being the only 
parameter; observations of a supernova in a metal poor region might show how 
the mass loss depends on metallicity. 

Subject headings: stars: circumstellar matter — stars: mass loss — supernovae 
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1. INTRODUCTION 

SNe IIP (Type IIP supernovae) may be relatively straightforward to interpret in terms 
of stellar evolution models. They are characterized by a long, plateau light curve that 
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can be attributed to the H envelope of the progenitor star. The progenitor stars have not 
undergone strong mass loss and are consistent with ~ 8 — 25 M stars that have evolved 
as single stars (e.g., Heger et al. 2003). Stellar evolution models, models of supernova light 
curves, and observations of supernova progenitors all suggest that these stars explode as red 
supergiants. 

Although only a fraction of the stellar mass is lost, red supergiants are known to have 
slow winds that provide a target for the expanding supernova. The properties of these winds, 
based on empirical fits to stellar observations, are an important input to stellar evolution 
models. The supernova interaction with the wind generates a hot, shocked region that can 
be observed by X-ray emission and by radio synchrotron emission from relativistic electrons 
(Weiler et al. 1986; Chevalier & Fransson 2003). SNe IIP provide an excellent place to 
study the interaction because asymmetries may not be significant and rapid evolution is 
not expected before the supernova. The main parameters determining the mass loss and 
explosion properties are expected to be the mass and metallicity of the progenitor star, 
although rotation could also play a role. 

The expected mass loss rates and hydrodynamic interaction are described in § 2. The 
resulting radio and X-ray emission properties are discussed in the context of observations 
of SNe IIP in § 3 and § 4, respectively. The relation of progenitor mass, deduced from the 
recent identifications of the stellar progenitors of SNe IIP, to mass loss is treated in § 5. 

2. HYDRODYNAMIC INTERACTION 

Calculations of massive star evolution typically include mass loss by using a parametrized 
fit to observed stellar mass loss rates. A number of evolutionary calculations (e.g., Schaller et 
al. 1992; Heger et al. 2003) use the parametrizations of de Jager, Nieuwenhuijzen, & van der 
Hucht (1988) or Nieuwenhuijzen & de Jager (1990). An older parametrization specifically for 
cool stars is that of Reimers (1977). Fig. 1 shows the values of M at the time of the explosion 
using the evolutionary tracks of Schaller et al. (1992), including the values given by Schaller 
et al. based on de Jager, Nieuwenhuijzen, & van der Hucht (1988) and values derived from 
the expression of Reimers (1977). The uncertainty in the mass loss rates is estimated by the 
authors to be about a factor of 2, which is comparable to the difference between the two 
parametrizations. More recent evolutionary tracks, like those given in Smartt et al. (2004), 
are in close agreement with those given by Schaller et al. (1992). Although the exact mass 
limits over which SNe IIP occur are not known, it can be seen that M varies over about an 
order of magnitude. 
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In these models, the progenitor star does not enter a 'superwind' phase at the end of its 
life. Eldridge & Tout (2004) find that the stellar luminosity can rise at the second dredge-up 
at the end of a star's life, presumably giving a rise in the mass loss rate, perhaps driven 
by thermal pulses. The upper initial mass limit for which the dredge-up occurs depends 
on convective overshoot: 8 M Q limit for overshooting and 10 M Q limit for no overshooting. 
If much of the envelope is lost, the supernova may not be of Type IIP. Chugai (1997) has 
discussed how supernovae with dense wind interaction may have progenitors near the low 
mass limit. 

The above results are for solar metallicity, Z m 0.02, where Z is the heavy element mass 
fraction. Evolutionary calculations typically assume that M oc Z a5 (Schaller et al. 1992; 
Heger et al. 2003), although this dependence does not have a sound basis in either theory or 
observations of cool stars. 

Another possible parameter for mass loss is stellar rotation. Meynet & Maeder (2000) 
find a significant dependence of the mass loss rates as a function of rotation. If we assume 
that most of the He-burning takes place in the red supergiant phase (as is true for the rotating 
models although the nonrotating low mass models spend some time as a blue supergiant) we 
find M = (0.84-1.6) x 10~ 6 M© yr" 1 for the 15 M models and (3.0-6.2) x 10~ 6 M© yr" 1 for 
the 20 M© models for non-rotating and 300 km s" 1 rotating models, respectively. Meynet 
& Maeder argue that the 300 km s -1 models are the more realistic. The values for the 
300 km s" 1 rotation models agree with the numbers in Fig. 1, while the non-rotating values 
are a factor of two lower. 

The important mass loss parameter for supernova interaction is the wind density, p w , 
determined by the parameter A, where p w = Ar~ 2 = M /Aitv w r 2 and v w is the wind velocity. 
The wind velocity for a red supergiant is typically 10 — 15 km s -1 . There is not a well- 
developed theory for the winds from red supergiants, but general considerations of winds 
typically have v w proportional to escape velocity, or v w oc M 1 / 2 //? 1 / 2 , where R is the stellar 
radius. For the range of cases shown in Fig. 1, the resulting variation in v w is < 10% about 
a mean value, with the 7 M© star having a high value and the 20 M© star having a low 
value. The variation in wind density is thus only slightly larger than that in M alone. 

For the structure of the supernova, we use the results of Matzner & McKee (1999) for 
the explosion of a red supergiant. These models assume that the outer envelope can be 
approximated as an n = 3/2 polytrope, which is accurate if the convection is efficient. The 
results for the exploded outer structure depend on q, the fraction of the total progenitor 
mass that is inside the H envelope. Since we are mainly concerned with stars that have 
not undergone extensive mass loss, we take q = 0.3. We have considered the interaction of 
the harmonic mean density profile of Matzner & McKee with a stellar wind typical of that 
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expected around a SN IIP and found that the density over the time of interest (~ 1 yr) is 
well approximated by the outer power law density region. The density in this region can be 
expressed as psn °c t~ 3 (r/t)~ n with n = 11.73. For this power law profile, the outer shock 
wave expands as R oc t( n ~ 3 )/( n ~ 2 ) oc t - 90 . If the wind interaction is considered in the thin 
shell approximation (Chevalier 1982), the shell velocity is given by 

\v wl J 51 VlOM ; \ 10 day,/ 1 ; 

where M_6 is the wind mass loss rate in units of 10~ 6 M Q yr -1 , v w \ is the wind velocity in 
units of 10 km s _1 , E§\ is the explosion energy in units of 10 51 ergs, and M e j is the ejecta 
mass. The temperature of the gas at the reverse shock front is T r = 1.2V^| keV, where K,4 is 
the shell velocity in units of 10 4 km s _1 . The model velocities can be compared to observed 
velocities in SNe IIP deduced from absorption lines, which give a minimum value for the free 
expansion velocity at the interaction point. Elmhamdi et al. (2003) find a maximum velocity 
of ~ 16, 000 km s _1 for the blue wing of H-alpha in SN 1999em at an age of 9 days. This 
velocity is difficult to determine accurately because of blending with other lines; however, 
the better defined absorption minimum is at 12,500 km s -1 , lending support to the above 
maximum velocity. 

Another consideration for the maximum ejecta velocity is that it may be limited by 
radiation losses at the time of shock breakout, independent of the circumstellar interaction. 
For a red supergiant with radius 500 R® and standard supernova parameters, equation (32) 
of Matzner & McKee (1999) gives a maximum velocity of ~ 13, 000 km s -1 , indicating that 
this could be a limiting factor. 

In the interaction model, the cooling time for gas at the reverse shock is equal to the 

age at an age 

t c w 20M_ 6 u- 1 1 Kl 5 - 2 days. (2) 
Substitution of equation (1) for V S 4 yields 

valid for a reverse shock temperature < 2 x 10 7 K, or for n = 11.73 and V s < 12, 000 km s' 1 
(Chevalier & Fransson 2003). Considering the dependence of M on mass (Fig. 1), it can 
be seen that stars with M > 20 M have an extended radiative period. However, most 
SNe IIP are expected to have progenitors of lower mass and to have shocks that are energy 
conserving except at very early times, before observations typically take place. The presence 
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of radiative cooling would be expected to give a boxy Ha profile at late times, as seen in 
SN 1993J (Filippenko, Matheson, & Barth 1994); such a profile has never been reported in 
a SN IIP. 



3. RADIO EMISSION 

Radio emission observed from supernovae can be attributed to synchrotron emission 
from the interaction region between the supernova and the surrounding wind (Chevalier 
1982). The mechanisms for particle acceleration and magnetic field amplification are not 
sufficiently well understood to predict the radio luminosity, L r . However, with the assump- 
tions that the particle and magnetic energy densities are a constant fraction of the total 
postshock energy density, we have L r oc (M / v w ) 1A ~ 2 for typical parameters, when the emis- 
sion is optically thin and at a fixed age (Chevalier 1982). This expression assumes similar 
supernova properties. 

At early times, the radio emission is absorbed; the most likely mechanisms for SNe 
IIP are FFA (free-free absorption) by ionized gas outside the forward shock front and SSA 
(synchrotron self-absorption). Optical depth unity to free- free absorption is attained at an 
age 

/ ■ \ 2/3 

where T cs ^ is the circumstellar temperature in units of 10 5 K. The value of T cs is uncertain. 
Calculations, mainly directed to Type IIL supernovae, by Lundqvist & Fransson (1988) find 
that for M_ 6 /v wl = 3 the temperature at the time of optical depth unity was ~3x 10 4 K, 
while for M_ 6 /v wl = 10 it was found that T cs ~ 1 x 10 5 K. These calculations used a simplified 
model for the shock breakout radiation, and may have underestimated the temperature 
somewhat. In the following we will use T cs = 1 x 10 5 K as a reference value, but we 
emphasize that this parameter is uncertain. 

A number of recent, well-observed SNe IIP are listed in Table 1. The distances to SN 
1999em and SN 2004dj are based on Cepheids; for SN 1999em, this distance is 50% larger 
than the distance determined by the EPM (expanding photosphere method) (Leonard et al. 
2003). The EPM distance to SN 1999gi is used; in this case, the Cepheid distances to nearby 
galaxies are only slightly larger (Leonard et al. 2002). The distance to SN 2002hh and SN 
2004et (both in NGC 6946) is based on an average of 3 methods (Li et al. 2005a). 

The positions of the SNe IIP, as well as other supernovae, in a peak radio luminosity - 
time of peak diagram are shown in Fig. 2, which is an update of Fig. 4 of Chevalier (1998). 
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The times of peak of the SNe IIP were determined from plausible fits to the available data. It 
can be seen that the SNe IIP cluster in a small region of the diagram. The constant velocity 
lines in the diagram are based on the assumptions that SSA is the dominant absorption 
mechanism and that there is energy equipartition between particles and fields (although 
the results depend only weakly on this assumption). If the velocity deduced in this way 
is compatible with the maximum velocities in the supernova, SSA is a plausible absorption 
mechanism. If the velocity is lower than the apparent maximum velocity, another mechanism, 
such as FFA, must be operating. However, the model assumes that p = 3 and that there is 
no cooling of radio emitting electrons; as shown below, these assumptions are probably not 
accurate. The position of the SNe IIP, somewhat below the 10, 000 km s" 1 line, is close to the 
velocities expected in SNe IIP, indicating that SSA may play a role but other mechanisms, 
such as FFA, may also be important. 

Depending on the parameters characterizing the supernova ejecta, the circumstellar 
medium, the magnetic field and the relativistic electrons, cooling of the relativistic electrons 
may be important. As we show below, this has important consequences for both the light 
curve and the spectrum. 

We first consider synchrotron cooling. Assuming that the magnetic field energy density 
is a fixed fraction e# of the the thermal energy behind the circumstellar shock, the ratio of 
the synchrotron time scale to the expansion time scale is given 




Therefore, in the restricted range M-q/v w i > 4 and e# > 0.1 synchrotron losses may be 
important. 

Inverse Compton cooling of the relativistic electrons by the strong photon flux from 
the ejecta can be important especially during the early phases. The ratio of the Compton 
cooling time scale and expansion time scale is 

tcomp ( U oX \ 1 / CB \ V4 / / y X -1/2 / t \ ^ 

t V2 x 10 42 erg s~V \0.l) \v wl J V10 GHz/ V 10 days J ' 

(6) 

where Lboi is the bolometric luminosity of the supernova. Depending on Lboi, M, and the 
highly uncertain value of e#, inverse Compton cooling may or may not be important. In 
general, low values of e B and M favor inverse Compton cooling. 

There are few well determined bolometric light curves of Type IIP supernovae. Obser- 
vations of the normal Type IIP SN 1992H show that L bol decreased from ~ 3 x 10 42 erg s~ x 
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to ~ 1.3 x 10 42 erg s^ 1 during the first ~ 50 days, and then stayed constant up to ~ 100 
days, when it dropped (Clocchiatti et al. 1996). There is, however, a considerable range in 
luminosities within the IIP class. Hamuy (2003) finds that the absolute magnitude during 
the plateau phase varies between My = —15.37 and My = —18.57, with My = —17.41 
for SN 1992H. From modeling of SNe IIP, Eastman et al. (1994) find that at day one 
Lboi ~ 2.5 x 10 42 erg s -1 , and during the plateau phase Lboi ~ (1 — 2) x 10 42 erg s -1 , 
in good agreement with SN 1992H. Chieffi et al. (2003) find a similar luminosity for SNe IIP 
during the first 50 days. The energy and envelope mass of the models are, however, adjusted 
to fit the observations, and are not calculated from first principles. Adopting a luminosity 
of L bol ~lx 10 42 erg s -1 , we find that for e B < 0.1 inverse Compton cooling is likely to be 
important during most of the plateau phase, except for very high mass loss rates or shock 
velocities. 

Because t comp oc e^ 4 and t sync h oc e B 3 ^ the Compton timescale is likely to be short 
in the cases where the synchrotron timescale is long, and vice versa. In Figure 3 we show 
how the cooling timescales as function of the parameter €bA* = €bM- 5 /v w i. Because both 
the synchrotron and Compton cooling timescales are proportional to (t/z/) ' 5 we plot the 
scaled cooling timescale t coo ;(tio/ ^io) ' 5 /t, where t\o = t/10 days and u w = u/10 GHz. The 
remaining parameter is then / 42 = (Lb o i/10 42 erg s _1 )(V s /15, 000 km s _1 )~ 2 . 

As an example, for M-$/v w \ — 1, 6b — 0.1, v — 10 GHz, Z42 = 1 and t = 10 days we 
get £ CO mp A ~ 1 an d ^synch/t ~ 0.3. For the same parameters, but with e# = 10~ 3 , we instead 
obtain t comp /t ps 0.3 and t sync hA ~ 10. In both cases the electrons are cooling, but in the 
former case by synchrotron cooling and in the latter by Compton cooling. 

The fact that cooling is likely to be important in many cases for the relativistic electrons 
means that the spectral index of the radio emission (flux oc v~ a ) is expected to be steep 
with a ~ 1.0, if the particles are injected with an energy spectral index p ps 2. As cooling 
becomes less important at late times, the spectrum should flatten toward a ~ 0.5. This is 
one of the important diagnostics of cooling. 

To show the possible range of radio light curves, we have calculated radio light curves 
for four different cases, two where cooling is important, one 'minimal cooling' case, and 
one purely adiabatic case. All calculations have p = 2.2, M- e /v w i = 5, n = 11.73 and 
V a = 15, 000 km s- 1 at 10 days. For the bolometric light curve we take that of SN 1992H 
and scale it by a factor / bo i. For each value of /boi we vary e# and e r to give a peak flux of 
2 mJy at 4.9 GHz at 5 Mpc, typical of the observed Type IIP light curves. Figure 3 and 
equations (5) and (6) show that one can have two qualitatively different cases of cooling. For 
high values of e r (and therefore small values of 6b) Compton cooling dominates. Conversely, 
large values of 6b imply synchrotron cooling. We therefore calculate two models typical of 
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these cases, one, model C, with / bo i — 1, 6b — 0.0013 and e r = 0.05, and one, model S, with 
/ bo i = 1, e B = 0.15 and e r = 0.001. 

To find parameters relevant for a realistic non-cooling case is non-trivial. To minimize 
inverse Compton losses, and still be within the range of reasonable bolometric luminosities, 
we take /boi = 0.05 and V s = 20, 000 km s -1 . During the plateau phase this corresponds to 
Lhoi ~ 1 x 10 41 erg s _1 , giving Z42 = 5.6 x 10~ 2 . From Figure 3 we see that minimum cooling 
then occurs for €bM_ 5 /v w i ~ 2 x 10~ 3 . With M_ 5 /v w i = 1, to reproduce the approximate 
turn-on time we require e r ~ 1 x 10~ 2 to give the approximate flux level. We refer to this as 
model M. Finally for comparison, we show one model, model A, where we have artificially 
put both synchrotron and inverse Compton cooling to zero. In Figure 4 we show the resulting 
light curves for 5.0, 8.46, and 22.5 GHz. 

The most interesting aspect of these light curves is the variation in the optically thin 
decline rate. While the flux in the adiabatic model, model A, falls off as F v oc £-(p+ 5 ~ 6m )/ 2 — 
t~ - 8 , the others have a flatter evolution. Here m — (n — 3)/(n — 2) = 0.9 for n = 12. This 
is especially true for Model C, which is dominated by Compton cooling of the electrons. 
Because the Compton cooling rate depends on the bolometric luminosity of the supernova, 
the radio light curve reflects the shape of the bolometric, primarily optical, light curve. In the 
limit that t comp <C t the stationary solution of the kinetic equation for the relativistic electron 
distribution (Fransson & Bjornsson 1998), together with the expression for the synchrotron 
flux, shows that 

F„ oc L hoX (t)-H m -v' 2 v- p ' 2 . (7) 

This shows that the synchrotron flux in this limit is inversely correlated with the bolometric 
luminosity. Further, for m ~ 1 and p m 2 one obtains a flat light curve in the plateau phase 
of the bolometric light curve. Note also the dip in the numerical light curves at ~ 100 days, 
seen especially at the highest frequencies. This coincides roughly with the time when the 
electrons become adiabatic, which at 5.0 GHz occurs at ~ 80 days for these parameters. 

In the synchrotron cooling dominated model S, the cooling has in contrast a smooth 
evolution, t s oc t 1 / 2 . Consequently, the radio light curve does not show any features like the 
Compton dominated case. Finally, we note that the minimum cooling case, model M, has a 
flatter decline than the purely adiabatic case, showing that cooling even for this case cannot 
be neglected. 

The only SNe IIP with sufficient radio data available to compare to our models are SN 
2004et and SN 2004dj. In other cases, there are some data on the radio rise, which can be 
used to estimate the circumstellar density from FFA. 
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3.1. SN 2004et 

The data on SN 2004et, shown in Figure 5, represent the most complete time coverage 
of any SN IIP. The data are from Stockdale et al. (2004c) before day 20, with additional data 
at 4.99 GHz after day 20 from Beswick et al. (2004) and Argo et al. (2005). Some additional 
upper limits at early epochs exist for the 4.99 GHz observations, but these do not add any 
strong constraints. Because of the good time coverage we have done a more detailed analysis 
of this supernova than for the other Type IIPs. 

The explosion date of SN 2004et is well constrained to be 22.0 September, 2004, within 
one day (Li et al. 2005a). The light curve of SN 2004et showed that it was a Type IIP 
supernova, although it showed some differences from typical Type IIP events (Li et al. 
2005a). For the simulations below we use the bolometric light curve of SN 1992H, scaled 
with /boi = 0.5. From the spectra of Li et al. (2005a) we find a maximum velocity of 
~ 14, 200 km s" 1 for the red wing of the Ha line at days 9 and 20. Because of occultation 
and other effects this is only a lower limit, and in the following we will use an expansion 
velocity of 15, 000 km s _1 at 10 days. The relative turn-on times between the different 
frequencies depend on the value of n, and we here take n — 10. A smaller value may give a 
better agreement with the observed turn-on times, but would also result in steeper decays 
in the optically thin parts of the light curves. All models have an injected electron spectrum 
with p = 2.2. Cooling considerably steepens the integrated spectrum. 

As discussed above, we vary the value of e#, and consequently e r , to cover both a 
Compton dominated and a synchrotron dominated case. In Figure 5 we show two sets of 
light curve fits with either large or small values of e# and e r , e# = 0.2 and e r = 0.001, and 
e B = 0.0013 and e r = 0.04, respectively. The value of M_ 6 /v wl is determined by the turn-on 
time, but, as pointed out before, also depends on the uncertain temperature of the CSM. In 
the first case we obtain M-q/v w i = 9 T c 3 / 5 4 , and in the second M-e/v w i = 10 T c 3 / 5 4 . 

The two cases result in fairly similar light curves, both strongly affected by cooling during 
the first ~ 100 days. The main difference is that for < 20 days the Compton dominated case 
shows a lower peak flux for the higher frequencies compared to the lower. This is caused by 
the decreasing Compton cooling rate with time, which affects the higher frequencies, which 
first become optically thin and also have a higher Lorentz factor of the electrons. In addition, 
as we have already found, at later epochs the light curves are flatter than for the synchrotron 
dominated case. The considerable error bars of the observations, especially at early time, 
make it difficult to favor one of these cases over the other. Because of the difference in 
decay slope it is of great interest to pursue these observations to as late times as possible 
to check these alternatives. The frequency dependent dip in the light curve at ~ 100 days 
in the Compton case would be an interesting signature. However, effects connected with a 
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gas density gradient different from a simple p oc r 2 wind may also influence the late time 
behavior and high observational accuracy would be required to detect the dip. 

In the e B = 0.0013, e r = 0.04 model, the optical depths to synchrotron self-absorption 
and free-free absorption are similar, while in the e B = 0.2, e r = 0.001 model the SSA optical 
depth is larger than 1 when it becomes thin to FFA. Therefore, especially in the latter 
model the light curve and spectrum are affected by SSA, giving a slower, more gradual rise 
compared to pure FFA. 

3.2. SN 2004dj 

A spectrum of SN 2004dj at the beginning of August 2004 indicated a normal Type IIP 
supernova (very similar to SN 1999em) with an age of about 3 weeks (Patat et al. 2004). We 
take 1 August, 2004 to correspond to an age of 21 days. The radio data in this case include 
a point at 8.4 GHz (Stockdale et al. 2004a), two at 1.4 GHz (Chandra & Ray 2004), and 
more extensive coverage at 4.99 GHz (Beswick et al. 2005). 

As in the case of SN 2004et, there is no single unique combination of e B and e r that 
fits the observations. In Fig. 6 we show two examples, one with small e B and large e r , 
and one for which the opposite is the case. The mass loss rate parameter is rather well 
constrained by the rising 1.4 GHz points, due to decreasing absorption. The observed sharp 
rise between days 33 and 43 is suggestive of FFA, although SSA cannot be ruled out. In 
fact, the e B = 0.01, e r = 0.0009 model is dominated by SSA, while FFA is dominant for the 
e B = 0.00075, e r = 0.01 model. We find that M_ 6 /v wl = 2-3 T c 3 / 5 4 for the two models. 

3.3. SN 2002hh 

Radio data on SN 2002hh are available from Stockdale et al. (2004b) and Chandra, Ray, 
& Bhatnagar (2003) for data at 1.4 GHz. The supernova was discovered on 31.1 October, 
2002 and a spectrum on 2 November, 2002 showed that it was a very young, highly reddened 
Type II supernova (Filippenko, Foley, & Swift 2002). We take the age at the time of discovery 
to be 2 days. There are no published optical light curve observations, as yet, that identify 
SN 2002hh as a plateau type. However, optical spectra of the supernova in the nebular phase 
show that it has slow moving H (Mattila et al. 2004); the mixing of H to low velocities is 
expected to require a massive H envelope compared to the core, as occurs in a SN IIP. Also, 
its radio properties are similar to other SNe IIP. 

The radio rise suggests optical depth unity to FFA at 1.4 GHz on day 62. Table 2 lists 
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the resulting value of M-qv~\T~^ 4 . This quantity is oc V^/ 2 ; the value V S 4 = 1.5 at 10 
days has been taken and allowed to evolve V S 4 oc t~ 01 . At the relatively low circumstellar 
density for SN 2002hh and other SNe IIP, Compton heating to ~ 5 x 10 4 — 10 5 K occurs and 
recombination is not important (see above, Lundqvist & Fransson 1988). 

A different picture for the circumstellar medium of SN 2002hh comes from its infrared 
detection which Barlow et al. (2005) find is due to dust in circumstellar mass loss ~5x 10 17 
cm from the supernova based on the infrared spectrum and luminosity. A similar radius for 
the dust can be obtained by considering that the duration of the infrared emission should be 
comparable to the light travel time across the emitting region and the total radiated energy 
should be approximately that radiated optically by the supernova. Barlow et al. (2005) 
estimate the dust mass to be 0.10 — 0.15 M , corresponding to a gas mass of 10 — 15 M , 
and attribute a significant part of the Ay fa 5 extinction toward the supernova to this dust. 
Supernova dust echoes were first inferred around SN 1979C and SN 1980K for which Dwek 
(1983) found circumstellar optical depths of 0.3 and 0.03 respectively. Scaling with mass 
loss density and using our value of M/v w leads to an optical depth < 0.01 for SN 2002hh, 
which would not contribute a significant extinction. An estimate of the mass loss rate from 
the mass and size of the dust emitting region also leads to a rate ~ 10 2 times larger than 
that estimated from the radio emission. 

In principle, there is no contradiction between these observations; the radio emission 
samples a region inside of the dust emitting region. There may be a rise in the radio emission 
in 15 — 20 years. However, the evidence for the loss of 10 — 15 M & of gas is difficult to reconcile 
with the slow moving H in the supernova, which would seem to require a massive H envelope 
at the time of the explosion. 

3.4. Other Supernovae 

SN 1999em was a well-observed, normal Type IIP supernova, including evidence for 
slow H in the nebular phase (Elmhamdi et al. 2003). The radio and X-ray data on it have 
already been discussed in terms of circumstellar interaction by Pooley et al. (2002). The 
peak in the light curve is suggested by the comparable fluxes at 1.4 and 4.9 GHz at days 
50 — 70 (see Fig. 2 of Pooley et al. 2002). If not for these, the data could only give an upper 
limit for the time of peak flux. However, the agreement of the results here with those from 
the other supernovae suggests that the flux peak was observed. An estimate of the mass loss 
rate is given in Table 2. 

The low luminosities of the detected SNe IIP show that the published upper limits that 
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have been set for radio supernovae (e.g., Weiler et al. 1989; Eck, Cowan, Branch 2002) are 
above the expected flux levels for supernovae of this type. However, an upper limit of 0.07 
mJy (8.46 GHz) for SN 2003gd on 17 June, 2003 (Stockdale 2003) is more stringent. The 
age of the SN IIP was ~ 87 days at that time (Van Dyk, Li, & Filippenko 2003). Assuming 
models similar to those discussed above, the limit is consistent with the expected flux of SN 
2004dj, but is a factor ~ 3 lower than the expected fluxes of SNe 2002hh and 2004et. These 
results suggest an upper limit M_ Q /v w i < 3. 



4. X-RAY EMISSION 

We first consider the thermal emission from the shocked region. The monochromatic 
luminosity of the circumstellar shock can be estimated from 

2 

^cs a n _ n r,35Am-0.08 7n-0.23 _-0.0116Ewir /To / \ ta-1 I f I I * 

dE 



6.9xl0 35 CT 9 -°- X?v 23e "°- 0116£kcv/T9 f ^ ) Vsti— ) 1 ( -fA—) 1 ergss^keV- 1 

V Vwi I \10 days J 



(Nymark, Fransson, & Kozma 2005). Here T 9 = T e /10 9 K, E^ e y is the photon energy in keV, 
/ is the relative thickness of the circumstellar shock region, and ( = [1 + 2n(He)/n(H)]/[l + 
4n(He)/n(H)] ps 0.86 for solar composition. For n = 7 we have / ~ 0.3, and for n > 12, 
/ 0.22 (Chevalier & Fransson 1994). The electron temperature, T e , is uncertain because 
of the unknown importance of collisionless heating. If only Coulomb collisions are important 
5 x 10 8 < T c < 1 x 10 9 K. The luminosity from equation (8) should be added to the flux 
from the reverse shock below. 

If the reverse shock is in the nonradiative regime and free-free emission dominates, which 
is a good approximation for T > 5 x 10 7 K, the X-ray luminosity at ~ 1 keV is given by 
(Fransson, Lundqvist, & Chevalier 1996) 

^ = 2.0xlO-C(n-3)(n-4)^«-e— /- (^)V (^)"' ergs S - 

(9) 

The main variation of X-ray luminosity with wind density is the (M_q/v w i) 2 term. Both 
equations (8) and (9) assume that the ingoing half of the X-ray emission is absorbed by the 
ejecta. In general, the emission from the reverse shock dominates in the keV region, unless 
the shock is radiative so that absorption by the cool dense shell is important. 

For T < 5 x 10 7 K, line emission becomes important and a more detailed calculation 
of the spectrum of the reverse shock is required. For this we have used the code described 
in Nymark, Fransson, & Kozma (2005) which calculates the ionization and emission for a 
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given electron temperature and density. Because the reverse shock, for the mass loss rates 
discussed here, is nonradiative, we assume that the ion temperature and density are given by 
the self-similar structure for the particular value of n. The electron temperature is calculated 
assuming Coulomb collisions only. For the reverse shock, this is enough to bring electrons 
and ions into equipartition. The ionization is calculated from the time a given element is 
shocked to the observed epoch, using the full time dependent equations. Because only a 
fraction of the hydrogen envelope has been lost for the Type IIP supernovae, we assume 
solar composition. 

In Fig. 7 we show the resulting X-ray spectra for n = 7 and 12, and V s = 10, 000 km s _1 
and V s = 15, 000 km s -1 . The mass loss rate is assumed to be 5 x 10~ 6 M Q yr _1 (for v w = 
10 km s -1 ). Because dL/dE oc (M_ 6 /v w i) 2 as long as cooling is not important, the spectra 
can be scaled to other mass loss rates, using this relation. The most interesting thing to note 
from Fig. 7 is the increasing importance of line emission as the reverse shock temperature 
decreases. For the models shown, T rev = 8.5 x 10 7 K for n = 7, V 8 = 10,000 km s _1 ; 
T rev = 1.9 x 10 8 K for n = 7, V s = 15,000 km s -1 ; T rcv = 1.7 x 10 7 K for n = 12, 
V s = 10,000 km s _1 ; and T rcv = 3.8 x 10 7 K for n = 12, V s = 15,000 km s^ 1 . While the 
line complex at ~ 1 keV is especially strong in the n = 12 models, it is absent in the n = 7 
models, where the only lines present are the Fe K lines at ~ 7 keV. The feature at ~ 1 keV 
is due to Fe XXI-XXIV, O VIII, and Ne IX-X, in this order. 

In Fig. 7 we have not included any absorption by the cooling gas behind the reverse 
shock, the circumstellar medium, or interstellar gas. The first two of these depend on the 
mass loss rate assumed. The total column density behind the reverse shock is 

*U»1.0xKP(„-4) ^^y^-'an-. (10) 

In the case that the gas is cooling (i.e. for M_6 > 5 and for early epochs, see eq. [3]) the 
emission from the reverse shock will be absorbed below ~ 1.2(iV rcv /10 22 cm~ 2 ) 3//8 keV. Even 
if the gas is not cooling and T rcv <2x 10 7 K, incompletely ionized atoms like Mg, Si, and 
Fe may in principle contribute to absorption above ~ 2 keV (e.g., Krolik & Kallman 1984). 
However, the opacity at 2 keV is only ~ 2.5 x 10~ 24 cm 2 for ~ 10 7 K, and absorption is 
unlikely to be important for reasonable parameters unless the gas is cooling. 

Both equation (9) and Fig. 7 show that the luminosity is sensitive to n. We found from 
the models for the radio emission for SNe IIP that n > 12 is indicated (§ 3). However, for 
the best observed X-ray SN IIP SN 1999em, Pooley et al. (2002) suggested n = 7 - 9 based 
on studies of the optical spectrum and the moderately high temperature (kT = 5 keV) found 
from the X-ray spectrum. The result was M_ 6 /v wl = 1 — 2; a higher value of n would lower 
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the mass loss density, bringing it further from agreement with the value deduced from the 
radio observations (Table 2). 

X-rays may also be produced by inverse Compton emission. As we have seen in § 3, 
inverse Compton cooling by the photospheric photons may be important. These photons 
are up-scattered to ~ 3.i7 2 (T e ff/10 4 K) eV, where 7 is the Lorentz factor of the relativistic 
electrons and T e g the effective temperature of the photospheric emission (Felten & Morrison 
1966). Using the expression for the inverse Compton emissivity in Felten & Morrison (1966), 
one can find the X-ray luminosity, which for the interesting case p = 3 (cooling important) 
takes a simple form, 

dLic QQ in38 p-i (M-6\„ ( L hol (t) \ ( t y 1 _ M x 

— « 8.8 x 10 e rlmm E kcV \ — j V s4 [ ^^^ ) ) ergs s keV . 

(11) 

Here 7m j n is the minimum Lorentz factor of the relativistic electrons. For other values of 
p, dLic/dE oc E^te -1 )/ 2 . The distinguishing property is the power law form, which is the 
same as that of the optically thin radio spectrum, and the correlation of the time evolution 
with that of the bolometric luminosity. It is unaffected by any absorption by cool gas 
behind the reverse shock. An observation of an inverse Compton component would allow 
a determination of the important parameter e r 7 min , and then from the radio observations 
€b- The stronger dependence of thermal emission on the circumstellar density implies that 
thermal emission is relatively more important at high M/v w and inverse Compton at low 
M/v w . At high mass loss rates, when reverse shock cooling becomes important, the inverse 
Compton component may again dominate, especially at early epochs. 

The observed X-ray luminosities are given in Table 1. The data are all from Chandra 
and refer to an age of ~ 30 days. SN 2004dj has an X-ray luminosity slightly higher that 
SN 1999em and a radio luminosity that is slightly lower (Table 1). This suggests that they 
have comparable wind densities. Both the radio turn-on and the X-ray luminosity of SN 
2002hh indicate that the circumstellar density is somewhat larger in this case. The low 
X-ray luminosity of SN 1999gi suggests a smaller density in this case (by a factor ~ 3), if 
the supernova properties are similar to the other cases. 

Because of the comparatively low fluxes, the spectral information from the observed 
supernovae is limited. A detailed comparison with the model spectra, or even a discrimina- 
tion between thermal and non-thermal spectra, is therefore not possible. Our calculations 
do, however, show that deeper observations of nearby SNe IIP can give very important 
constraints on both the mass loss and the shock physics. 
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5. PROGENITOR MASS AND MASS LOSS 

Progenitor mass estimates are from observations of the region of the supernova before 
the supernova occurred. Studies of this kind have been possible for a number of nearby SNe 
IIP (see Table 1). In the cases of SN 2003gd and SN 2004et (Li et al. 2005a), progenitor 
stars have been tentatively identified and, in the case of SN 2004dj, a mass was determined 
from the properties of a tight cluster at the position of the supernova (Maiz-Apellaniz et al. 
2005). Both SN 2004dj and SN 2004et are estimated to have progenitor masses ~ 15 M for 
which standard mass loss prescriptions give M_ 6 = 1.5 — 3; the range is increased to 1 — 7 
when the mass uncertainty for SN2004et is included. The mass loss rate deduced from radio 
observations of SN 2004dj is consistent with expectations, while that for SN 2004et suggests 
a mass ~ 20 M . However, for SN 2004et Li et al. (2005a) find that their tentative progenitor 
is a yellow, rather than red, supergiant and may have some similarity to the peculiar SN 
1987A and SN 1993J, although less extreme. These supernovae had very different radio 
properties from the normal SNe IIP (see Fig. 2), which have properties close to those of SN 
2004et. 

If the progenitor mass of SN 2003gd was 8 — 9 M , the mass loss density was a factor of 
~ 3 lower than for a 15 M & star (Fig. 1), leading to a radio luminosity 5 — 10 times smaller 
at the same age. This is consistent with the upper limit on the radio luminosity and shows 
the difficulty of observing low mass SNe IIP at radio wavelengths. The observations must 
be undertaken at an early time and with high sensitivity. The recent Type IIP SN 2005cs in 
M51 is also estimated to have a relatively low mass progenitor, 8 ± 1 M (Li et al. 2005b) or 
9^2 M Q (Maund, Smartt, & Danziger 2005). It has not been reported as an X-ray or radio 
source. 

The initial mass deduced for SN 1999em by Smartt et al. (2003) of < 15 M is lower 
than that found for SN 2004dj by Maiz-Apellaniz et al. (2005), but Leonard et al. (2003) 
found a mass < 20 M for SN 1999em. The value of M_ 6 /v w i (Table 2) is comparable to 
that for the other SNe IIP and is consistent with stellar evolution and typical red supergiant 
mass loss rates for a 15 M initial mass star. 

SN 1999gi was considerably less luminous in X-rays than SN 1999em (Table 1), implying 
a lower value of M_ 6 /v w i by a factor ~ 3 for comparable supernova properties. Considering 
the relation between M and M shown in Fig. 1, a drop in M by a factor 3 corresponds to 
going from a 15 M initial mass star to an 8 M & initial mass star. This is consistent with 
the limit < 12 M deduced for the progenitor of SN 1999gi. 

An independent way to estimate the mass of an exploded star is from the properties of 
the supernova. From analyses of the light curves of the supernovae, ejecta masses for SN 
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1999em have been estimated as 15.0 M (Nadyozhin 2003, for d = 11.1 Mpc), 27±§ 4 M 
(Hamuy 2003, for d = 10.7 Mpc) and, for SN 1999gi, 18.7 M (Nadyozhin 2003, for d = 11.8 
Mpc), 43±?4 M (Hamuy 2003, for d = 9.0 Mpc). With this method, the mass estimate 
decreases as the distance estimate increases. To obtain the initial mass, the mass of any 
compact remnant and mass lost during the evolution must be added to the ejecta mass. It 
can be seen that the masses obtained in this way are larger than the masses obtained from 
the mass loss or from the direct progenitor observations, although Nadyozhin's value for SN 
1999em is fairly close to the value deduced by these other methods. The method is, however, 
sensitive to several factors. In particular, hydrodynamic mixing, as was found to be very 
important for the light curve of SN 1987A, may change the mass estimates considerably (e.g., 
Chieffi et al. 2003). More complex supernova models, including models of the spectra, are 
needed to see whether better agreement can be found between these methods. The current 
work makes specific predictions for SN 1999em and SN 1999gi, which have only upper limits 
from the direct progenitor observations. 

Table 1 includes estimates of the metallicity of the regions in which the supernovae 
occurred. The metallicities are either solar or within a factor of 2 of solar; this variation is 
expected to have a minor effect on the present discussion if the estimated Z 1 / 2 dependence 
of mass loss is correct. If SNe IIP are discovered in highly metal poor regions, there is the 
possibility of investigating the metallicity dependence of mass loss. 

As the sample of supernovae grows, it will be possible to compare the distribution of 
observed progenitor masses with that expected for a standard initial mass function. If we 
consider two mass groups, 8 — 13 M and 13 — 18 M , 4 of the 7 supernovae mentioned here 
(SNe 1999em, 2002hh, 2004dj, and 2004et) belong to the high mass group and 3 (SNe 1999gi, 
2003gd, and SN 2005cs) belong to the low mass group. For a Salpeter mass function, there 
would be a factor 1.7 fewer supernovae in the high mass group compared to the low mass 
group. At least part of the reason for the observed relative paucity of low mass progenitors 
may be that they give rise to subluminous supernovae. SN 2005cs had a peak absolute V 
magnitude of —15.6 (Li et al. 2005b), which is subluminous (see § 3). Also, it is possible 
that not all stars in the 8 — 13 M range end their lives as SNe IIP. The lower mass cut 
for supernovae is not well-determined, and stars in this mass range may end their lives as a 
type of supernova other than a IIP. 

In summary, data on the wind interaction and progenitor stars of a small number of SNe 
IIP are consistent with mass loss estimates and tracks in evolutionary models of ~ 8 — 20 M & 
single stars. The radio and X-ray emission from SNe IIP is relatively faint, but this emission 
may provide a check on the mass of the supernova progenitor. In addition, the fact that the 
radio emitting electrons are likely to be subject to energy loss affects the radio light curves 
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and gives constraints on the physical properties of the shock interaction. X-ray emission gives 
additional diagnostics and will be most useful if well-determined spectra can be obtained. 

We are grateful to Megan Argo, Rob Beswick, Chris Stockdale, Kurt Weiler, and their 
collaborators for making radio data available in advance of full publication. This research 
was supported in part by NSF grant AST-0307366, Chandra grant TM4-5003X, and the 
Swedish Research Council and National Space Board. 
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Table 1. Type IIP supernovae 



SN 


Distance 


Metallicity 


L x (0.5-8 keV) 


Progenitor 


Refs. 




(Mpc) 


(Z/Z e ) 


(ergs s" 1 ) 


mass (M Q ) 




1999em 


11.7 ±1.0 


1-2 


9 x 10 37 


< 15 


1,2,3 


1999gi 


11.1 


~ 2 


1.6 x 10 37 


< 12 


2,4,5 


2002hh 


5.5 ±1.0 




4 x 10 38 




6,7 


2003gd 


9 


~ 0.5 




8-9 


8,9 


2004dj 


3.2 


1 


1.5 x 10 38 


15, ~ 12 


10,11,12 


2004et 


5.5 ±1.0 






15^ 


6 



References: (1) Leonard et al. 2003; (2) Smartt et al. 2003; (3) Pooley et al. 2002; (4) 
Leonard et al. 2002; (5) Schlegel 2001; (6) Li et al. 2005; (7) Pooley & Lewin 2002; (8) Van 
Dyk et al. 2003; (9) Smartt et al. 2004; (10) Maiz-Appelaniz et al. 2004; (11) Wang et al. 
2005; (12) Pooley & Lewin 2004. 



Table 2. Results for radio supernovae 



SN 


t( Tff = 1) at 1.4 GHz 
(day) 




1999em 


52 


5 


2002hh 


62 


7 


2004dj 




2-3 


2004et 




9-10 
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j_ Q-7 I I I I I I I LJ 

10 

M (M ) 

Fig. 1. — Mass loss rates from red supergiants at the end of their lives as a function of 
initial mass. The stellar endpoints are from M — 7,9, 12, 15, 20, 25 M evolutionary tracks 
by Schaller et al. (1992), and the mass loss rates are based on de Jager et al (1988) (solid) 
and Reimers (1977) (dashed). 
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Fig. 2. — Peak radio luminosity vs. time of peak for radio supernovae. The 4 Type IIP 
supernovae (crosses) are near the middle of the plot. 
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Fig. 3. — Ratio of cooling time scale to adiabatic time scale for synchrotron cooling (single 
solid line) and for inverse Compton cooling (4 parallel lines). For Compton cooling each 
curve is labeled with the value of Z 42 = (L bol /10 42 erg s -1 )(V^/15, 000 km s _1 )~ 2 . 
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Fig. 4. — Radio light curves of four models with different cooling properties. The dotted line 
corresponds to a purely adiabatic model, the short dashed to a 'minimal cooling' model, the 
long-dashed curves to a model dominated by synchrotron cooling, and the solid curves to a 
model dominated by Compton cooling. See text for the parameters of each model. 
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Fig. 5. — Observed radio light curves for SN 2004et compared to models with two different 
combinations of €b and e r . Note especially the effects of Compton cooling for the upper 
panel, resulting in nearly flat optically thin light curves, as well as the dip at ~ 100 days 
seen at the high frequencies. 
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Fig. 6. — Observed radio light curves for SN 2004dj compared to models with two different 
combinations of e# and e r . 
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Fig. 7. — Calculated X-ray spectra for n = 7 and 12, and V s = 10, 000 km s _1 and V s = 
15,000 km s _1 . The mass loss rate is 5 x 10 -6 M Q yr _1 (for v w = 10 km s _1 ). Note the 
increasing importance of line emission as the reverse shock temperature decreases. For the 
highest temperature (n — 7 and V s = 15, 000 km s" 1 ) the only important lines are the Fe K 
lines at ~ 7.0 keV. 



